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Abstract 

The relativistic Maxwell-Boltzmann distribution for the system of A'^ events 
with motion in space-time parametrized by an invariant "historical time" r 
is considered without the simplifying approximation m? = M^, where M is 
a given intrinsic property of the events. The relativistic mass distribution is 
obtained and the average values of m and rm? are calculated. The average value 
of the energy in nonrelativistic limit gives a correction of the order of 10% to 
the Dulong-Petit law. Expressions for the pressure and the density of events 
are obtained and the ideal gas law is recovered. 
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1 Introduction 

This article is the natural continuation of the work of Horwitz, Shashoua and Schieve 
which in the present text will be cited as I. In that work a derivation of a 
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manifestly covariant relativistic Boltzmann equation for the system of events was 
made in the framework of a manifestly covariant classical and quantum mechanics 0. 
These events, considered as the fundamental dynamical objects of the theory, move 
in an 8A^-dimensional phase space. Their motion is parametrized by a continuous 
Poincare invariant parameter r called the "historical time" . The collection of events 
(called "concatenation" 0) along each world line corresponds to a particle in the 
usual sense; e.g., the Maxwell conserved current is an integral over the history of the 
charged event @] . Hence the evolution of the state of the iV-event system describes 
the history in space and time of an A^-particle system. The relativistic Boltzmann 
equation was used to prove the if-theorem for evolution in r. In the equilibrium limit 
a covariant form of the Maxwell-Boltzmann distribution was obtained. Since this 
distribution is the distribution of the 4-momenta of the events, = — = —p'^p^ 
is a random variable in a relativistic ensemble. In order to obtain a simple analytic 
result the authors restricted themselves to a narrow mass shell p"^ = —m? = — 
where M is a given fixed parameter, with the dimension of mass (an intrinsic property 
of the events), assumed to be the same for all the events of the system. The results 
obtained in this approximation are in agreement with the well-known results of Synge 
1^ from an on-mass-shell relativistic kinetic theory. 

As Hakim has remarked, some problems of relativistic statistical theory re- 
quire consideration of relativistic particles endowed with variable masses. Thus, in 
statistical cosmology the universe is constituted of a gas of particles with equal 
masses. The particles are, in turn, considered as being galaxies and possibly clusters 
of galaxies or stars. However, galaxies and clusters of galaxie^ do not appear to 
have identical masses. Therefore, a consistent approach to statistical cosmology as 
well as the statistical treatment of clusters of galaxies themselves should involve a 
mass distribution. 

Moreover, in the statistical bootstrap model of Hagedorn and Frautschi |^ 
for multiple production of particles in high energy reactions a mass spectrum of 
the asymptotic form p{m) ~ cm"exp(6m) (where a,b,c are constants) lies in the 
foundations of their theory and gives a good qualitative agreement with relevant 
experiments in high-energy physics. Miller and Suhonen have discussed a possible 
correlation of the grand canonical distribution function of ref.iin], characterized by 



the mass fluctuations, with the hadronic spectrum of the Hagedorn-Frautschi form. 

In this paper, we study such a relativistic system without the simplifying approx- 
imation = M^, i.e., we consider a relativistically covariant Maxwell-Boltzmann 
distribution with mass parameter m taking any value within the range < m < oo. 
We obtain the relativistic mass distribution by integration of the Maxwell- Boltz- 
mann distribution over angular and hyperbolic angular variables. Calculation of 
average value of the energy and study of its nonrelativistic approximation gives a 



number of clusters, such as Coma or Virgo, contains large numbers of galaxies of unequal 
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relativistic correction of the order of 10% to the Dulong-Petit law for the free-particle 
gas assuming that the mass distribution remains valid in this approximation of our 
theory. We remark that, as is well known |]TT|, the structure of the Galilean group, 
the symmetry of non-relativistic systems, implies that the mass of a particle must 
be a constant intrinsic property. We recognize, however, that the applicability of 
the Galilean group is an idealization of a world which seems to be more correctly 
described by the Poincare group. The result that we have found follows from equi- 
librium thermodynamics without imposing the geometrical restriction of the precise 
Galilean group to an infinitely sharp mass shell. 

By examining the energy-momentum tensor we obtain expressions for the pressure 
and the density of events and recover the ideal gas law previously obtained in I. 

Since there was no fundamental theory available to Hakim [0, his development 
of the theory followed the phenomenological methods initiated by Synge [5] . In this 
way he obtained a relativistic statistical mechanics containing many useful results 
|I3| . In particular, he used the Jiittner- Synge form (coinciding with the sharp-mass 
or low-temperature limit of the theory which we discuss here) to obtain a mass dis- 
tribution [6]. As we shall see, this construction leads to a normalization of the dis- 
tribution function, which differs from ours. The theoretical framework which we use, 
permitting a derivation of the relativistic Boltzmann equation from first principles, 
provides a condition on the equilibrium ensemble which leads to a well-defined mass 
distribution (independent of the Jiittner-Synge result, although consistent with it), 
and a normalization condition consistent with a quantum-mechanical positive definite 
density. 



2 Relativistic mass distribution 

We begin with the Maxwell-Boltzmann distribution function which we write down in 
the form (differing from one obtained in [1,(3.14)] by the signs in the exponent^]; we 

^These signs do not influence the requirements for the MaxweU-Bohzmann distribution because 
the logarithm of the distribution function 

ln/o(g,p) - A(p + pe)'+lnC(g) 

possesses the property [1,(3.10)] 

ln/o(g,p) = Xii<l,P) +X2{q,p) + ■■; 

where the quantities 

Xi{q,pi) + Xi('i,P2) 

are conserved in colhsions. For the sharp mass approximation = the difference between the 
two expressions manifests itseff only in the normalization factor. 
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also use the metric g^^ = (— , +, +, +) and q = q^, p = p^) 

/o(g,p) = C(g)e^(^)(^+^=)^A(g)>0. (1) 
The function (1) must be normahzed, according to [1,(2.14)], as 

niq) = J fo{q,p)d'p = C{q) J d'pexp{Aip + p,f}, (2) 

where n{q) is the total number of events per unit space-time volume in the system in 
the neighborhood of the point q. By introducing hyperbolic variables [1,(3. 16), (3. 17)] 

fi^ : m > 0, < ^ < TT, < < 27r, -oo < /3 < cx), 

we can rewrite (2) as follows: 

n(q) = C(q)e'^"'' f sinh^ f3 sin 6 dmd (3 dO dipe-^"''-^'^"'^'^^°'^ . (3) 

After some calculations [Q we obtain the normalization relation 



n{q) = C{q)e~^'-' ^^^{2,2; Ami). (4) 

Here \l/(a, b; z) is the confluent hypergeometrical function of z (ref. [15], p. 257, sec. 6. 6). 
After integration over angular and hyperbolic angular variables we obtain from (3) 
the expression 

n{q) = C{q) — / dmm^e"^'^ K^{2Amcm), (5) 

AttIc Jo 

from which we identify the mass distribution function 

/(m) = Dm^e-'^"''Ki{2Am^m), (6) 

where 

= {2,2; Am% (7) 
so that f{m) is normalized according to 

f{m)dm = 1. (8) 



In (5) and (6) Ki is the Bessel function of the third kind (imaginary argument), 
where, in general, 

Kiz) = ^e-/^Hi'\zz). 



With the help of this distribution one can obtain the local average value of an 
arbitrary function of mass (j){m): 

roo 

{(f){m))g = / (j){m)f{m)dm. (9) 
Jo 

Let us obtain, for example, the local average values of mass and mass squared in 
relativistic gas which represent the first two moments of the distribution (6): 

Stt v&(|,2;Amg) 



^ 2 vl/(3,2;Amg, 



More generally. 



Now, as in /, we define absolute temperature through the relation [1,(3.26)] 

2Am^ = — -, 13 

which implies that in thermal equilibrium AttIc is independent of q. Hence 

Stt r — ' 



(m) = "-^prn^ksT - (14) 



In the limit T — > it follows from the asymptotic formula for z ^ oo [jT5[ 

a(a + 1) ■ ■ ■ (a + n - 1)(1 + a - 6)(2 + a - 6) ■ ■ ■ (n + a - 6) - 



n=l 



vI/(a,6;^)~^-'^[l+^(-l)^ 
that 



(15) 

(m) - ^fc^T. (16) 



One can also obtain in this limit 

(m^) - S{kBTf. (17) 
Let us now calculate the first two moments of the distribution (1): 

= J ^ ^ CIS) 
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If we introduce the "free energy" of the relativistic gas through the relation 



I d'^pe^^P+P'^^" = e 
we will obtain the following expressions: 



AF 



(F' - 1)2 - 



A 



+ 9 



F'-l 
2A ' 



F' = 



dF 

drrii 



(19) 

(20) 

(21) 
(22) 
(23) 



The value of F, obtained from (20) and (4), is 

1 



Hi 



^{2,2- Ami)]. 



A simple calculation with the help of the relation (ref.[15], p. 257, sec. 6. 6) 

d 



gives 



— *(a, b] z) = -a*(a + 1, 6 + 1; 2;) 
dz 



iP )<i ''Pc^{2,2;Amiy 



^(4,4;^m; 



-p'cPc + 



gi"" *(3, 3; Ami) 



(24) 



(25) 



'*(2,2;Am2)"''^"'" ' A *(2,2;Am2)' 

As in I, we make a Lorentz transformation to the local average motion rest frame 
moving with the relative velocity 



u 



in order to obtain the local energy density. 
The rest frame energy is 



{E') 



{E)g p 



so that 



^^^^-'^^¥(2:2^ 



(26) 
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Using the asymptotic formula (15) and expression (16), we have for T —>■ 0, 

{E') - (m) = (4 - ^)kBT = ^hsT, (27) 

16 — Stt , , 

^^—^^1.1. (28) 

This result differs from the nonrelativistic result ^ksT. The coefficient 7 represents a 
relativistic correction, determined by the relativistic mass distribution, to the classical 
value |. The existence of such corrections was first shown by Horwitz, Schieve and 
Piron [10] in their work on study of relativistic Gibbs ensembles even in the case of 
small fluctuations of mass over its sharp value Mq. 

In the limit T — > 00 we use the asymptotic formulas for 2; — > (ref.[15], p. 262, 
sec. 6. 8) 

^{a,b;z) ^ ^^^~ ^K '-' + 0{\ z f''-^), Reb>2, b ^ 2, (29) 
i [a) 

<i,{a,b;z)^^^^^f^z'-' + 0{\]nz\), 6 = 2, (30) 
r(a) 

and obtain 

{E') ^ 2kBT, (31) 

the result obtained also in I. 

In this limit one can find from (10), (11) 



^ hmickBT 

W = y — 2 — ' ^ 

(m^) ^ 2mckBT. (33) 

To obtain the pressure and density distributions in our ensemble, as in I, we study 
the particle energy-momentum tensor defined by the density 

T'''iQ) = i:Jdr^S\q-q,{r)). (34) 

i 

Using the result of the previous work [1,(3.37)] 

{T'%q)), = TAvJd'pfo{q,p)^ (35) 



and the expression (25) for {p'^p^)q, we obtain 



TAvn{q) 
M 



^^(3^3^Am|) ^ , ^(4,4;Amg) 
A ^(2,2;Am2)^ ^"^'^ *(2, 2; Am2) 



(36) 
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In this expression T^v is the average passage interval in r for the events which pass 
through the small four-volume AV over the point q of R^. 

The formula for the stress-energy tensor of a perfect fluid has the form [1,(3.39)] 

(T-(,)>, = P9--(p + p)|;|^. (37) 

where p is the pressure and p is the density of energy at q. 
According to (24), 

, <^>' = ^, (38) 

hence 

^ T^ynjq) ^{3,3- Ami) .^q^ 
^ AM *(2,2;>lm2) 

and 

6TAyn(g)m^vI/(4,4;Amg) 
^ + M n^,2;Amiy 

To interpret these results, as in I, we should calculate the average (conserved) particle 
four-current, which has the microscopic form 

n<l)-i:J^^'i^-^ii^))dr. (41) 

i 

Using the result of / [(3.45), (3.59)] 

{J'{q)), = TAvJd'p^fo{q,p) (42) 

and expression (24) for (p'^)q, we obtain 

^ 2T^vn(q) M3,3;Aml) 
W WJ;. ^ ^'^*(2,2;Am2)- ^^"^^ 

In the local rest frame p(f = {rric, 0), 

0/ XX _ 2rAyn(g)mc^(3,3; Amg) 

^ ^(2,2;Am2)- ^^^^ 

Defining the density of paricles per unit space volume [1,(3.48)] as 

A^o(g) = (J°(?))„ (45) 
we obtain the ideal gas law [1,(3.49)] 



Now we shall show how the general expression for the distribution function (6) 
can be simplified in two limiting cases T ^ and T ^ oo. 

1. T — i> 0. In this case we can see from (17) and (27) that (m^) ~ [ksTY ^ 
{E')mc ~ mc{kBT). We therefore neglect = —p"^ in comparison to 2Appc in the 
exponent of the Maxwell-Boltzmann distribution (1); hence, we begin with 

ft>'\p,q) = Ciq)e-^"''^e'^^^^, (47) 
and after integration and normalization obtain 

f{m) = ^^^y^m^Ki{2Am,m). (48) 

It is easy to see that the distribution function (48) produces results for (m) and (m^) 
in this limit in agreement with (16) and (17): 

(m)° = ^knT, {my = 8(A:^T)^ 

More generally, 

{mr = r(^ + l)r(^ + 2){Am^)-' = T{^- + l)r(^ + 2){2ksTY, (49) 

which coincides with the low-temperature limit of (12). Similarly, if we introduce the 
"free energy" for the distribution (47), we obtain 

(p/^)0 = P^AkBT, (50) 
rric 

which coinsides with the low-temperature limit of (24). 

2. T —>■ oo. It is seen from (31) and (33) that (m^) and {E')mc are of the same order 
in T. But the argument of the function Ki in the expression of the mass distribution 

function (6) is of the lower order of T in this limit [(32)], i.e., 2Amcm ~ mi {kBT)2 . 
Since we use the asymptotic formula [1^ 



2 ' ' \2, 

which gives Ki{2Amcm) = l/2Amcm, and obtain the (normalized) distribution 

/°°(m) = 2^7716-^™'. (51) 

One can easy check that the distribution (51) gives results for (m) and (m^) in this 
limit in agreement with (32) and (33): 
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More generally, 



{m'r = r(- + = r(- + l)i2m,kBT)-^, (52) 



which coincides with the high-temperature limit of (12). 
In this way one can also obtain 



(pM)oo ^ £j^2kBT, (53) 



rrir 



the high-temperature limit of (24). 

Summarizing, we can write down the expression for the mass distribution function, 
taking into account the limiting cases considered: 



/(m) 



' {{2Amcf/2}m^Ki{2Amcm), T ^ 

{AA/mc'^{2,2] Aml)}e~'^"^^m^Ki{2Amcm), intermediate case (54) 
2Ame-'^"'\ T oo 



Synge [5], following phenomenological methods in his study of an on-mass-shell 
relativistic ensemble, used the normalization relation (the only normalization relation 
available to him) for the equilibrium distribution function for such a system^, 

fo{q,p)=G{q)e'^PP% (55) 

called the "Jiittner-Synge distribution function", 

/ fo{q,p)d'p = G{q) J dW"""^ = No{q), (56) 

where Noi^q) = {Jo{q))q is the density of particles per unit space-volume (45), the 
0-component of the average (conserved) particle four-current; in contrast to the nor- 
malization relation (56), we have the condition (2), 

fo{q,p)d^P = n{q), 

where n{q) is the density of events per unit space-time volume, a quantum-mechanical 
positive-definite density, which follows from the theory discussed in I. 
In the theory of Synge 

J9%4 Giq) 9$ 



{nq)),= ffo{q,p)-d'p 
J m 



m d{2Ap: 



his notation, 2Api^ — 
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where m is a (constant) mass of the particles of the system and (ref.[5], p. 35) 



$ = / e^^-d^p = 4^^^i(2-4r»cm)_ ^^^^ 

In this way he obtained the normahzation relation (ref.[5], p. 36) 

A'KGw?K2{2Amcm) 
^° 2A^^ ' ^^^^ 

differing from the one obtained in / [1,(3.18)] in the mass- sharp- approximation = 
for the same distribution function (55) within the framework of the theory dis- 
cussed in I: 

where G = C{q)exp{A{M + Mc)^} and Am is the fluctuation of mass around its 
sharp value M. 

Hakim introduced variation of mass into the Jiittner-Synge distribution and, using 
the normalization condition (58), found the mass distribution [6,(3.3)] 

f(m) = '^i^^^m''K2{2Amcm). (60) 
This distribution differs from our low-temperature limit (48) 

because it is obtained by assuming a different normalization of the initial Jiittner- 
Synge distribution. This is why the results of Hakim differ quantitatively from the 
results obtained within the framework of the theory discussed in the present paper 
in the low-temperature limit. 



3 Concluding remarks 

We have considered an equilibrium relativistic ensemble, described by the equilibrium 

relativistic Maxwell-Boltzmann distribution, with variable mass. Wc have found that 
the equilibrium state of such a system is characterized by a well-defined mass distri- 
bution (consistent with the Jiittner-Synge one in the sharp-mass or low-temperature 
limit), following directly from the Maxwell-Boltzmann distribution, by integration 
over angular and hyperbolic angular variables and satisfying a normalization condi- 
tion consistent with a quantum-mechanical positive definite density. 

The results of Hakim, who introduced variation of mass into the equilibrium rela- 
tivistic Jiittner-Synge distribution and found a relativistic mass distribution [6] , differ 



11 



quantitatively from the results obtained within the framework of the theory discussed 
in the present paper in the low-temperature limit (e.g., he did not obtain a relativis- 
tic correction to the Dulong- Petit law) because the initial Jiittner-Synge distribution, 
from which the mass distribution of Hakim is obtained, has a normalization differing 
from ours, based on the counting of world lines (particles) rather than events. 

The relativistic mass distribution obtained in this paper may find applications 
in relativistic cosmology and astrophysics (when dealing with clusters of galaxies or 
stars), in the statistical mechanics of hadrons, in physics of relativistic plasmas and 
relativistic liquids. 

In the present paper as well as in the previous one (/) the simplest case of relativis- 
tic equilibrium is treated. The nonequilibrium case may be treated separately, e.g., 
by application of the Grad method of moments adopted for the relativistic theory. 

The case of indistinguishable particles will be also treated separately fl^, in order 
to obtain their equilibrium relativistic distribution function and the corresponding 
equilibrium relativistic mass distribution. 
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